Derivative Expansion of One-Loop Effective Energy of Stiff Membranes with Tension 
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With help of a derivative expansion, the one-loop corrections to the energy functional of a nearly flat, stiff 
membrane with tension due to thermal fluctuations are calculated in the Monge parametrization. Contrary 
to previous studies, an arbitrary tilt of the surface is allowed to exhibit the nontrivial relations between the 
different, highly nonlinear terms accompanying the ultraviolet divergences. These terms are shown to have 
precisely the same form as those in the original energy functional, as necessary for renormalizability. Also 
infrared divergences arise. These, however, are shown to cancel in a nontrivial way. 
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CIh' Due to their small surface tension, fluid membranes (see Ref. [|l|] for reviews) are subject to strong thermal undulations. The 
energy of such a membrane is usually modeled by the local expression [HsJ 
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I where dS are the surface elements, while H and K denote (twice) the mean and the Gaussian curvature of the membrane surface, 
T-j. . respectively. In terms of the principal radii i?i and R2 of curvature these are 7J = 1/Ri + I/R2 and K = l/i?i_R2- The param- 
(yQ eter /ig in ([l|) is the surface tension, kq the bending rigidity, and Rq its Gaussian counterpart. The geometric quantities appearing 
7—i in the integral (|l]) are invariant under translations and rotations in space. They are also independent of the parametrization of the 
surface. We ignore a possible spontaneous curvature term linear in H. The energy of a physical membrane contains all higher 
CD powers in the principal radii of curvature, but these are iiTelevant at large length scales. (In the language of renormalization 
group analysis, the first term is relevant, the second and third are marginal.) 

The statistical behavior of fluctuating membranes was first studied by Helfrich ^ using only curvature terms. We allow for an 
additional surface term in (|l]) because flucutations arising from the curvature terms not only renormalize the bending rigidities, 
but also the tension. The thermal fluctuations soften the bending rigidity at large length scales, reducing it from the bare value 
Ko as follows 



Ch , ftoff = Ko- -7-3 In(AL), (2) 



where j3 is the inverse temperature, A is an ultraviolet momentum cutoff of the order of the inverse microscopic length scale a 
. !^ given by the length of the molecules, whereas L is an infrared cutoff determined by the finite size of the membrane. Various 
authors derived different values for a, first a — 1 [Q-^ was obtained, later a = 3 [[7[-p^. The second result has also been found 
■ in computer simulations jll]]. For either value of a, the rigidity disappears at length scales larger than the persistence length [ p^ 

^ ~ aexp f — /3koV (3) 
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beyond which the normal vectors of the surface become uncorrelated — the surface looks crumpled. More recent calculations 
1^ ] suggest the value a — —1, implying a stiffening instead of a softening of the bending rigidity. This new result was argued 
to arise from the use of another integration measure which respects the incompressible-fluid nature of the membrane from the 
outset. This is in contrast to previous studies of in-plane fluid and elastic effects [ [l4} | which did not show any change in the 
value a — 3 (they only enter at the two-loop level [|l5p). 

The renormalization of the Gaussian rigidity Rq was first calculated in to have the same form as in (||), but with a ^ a = 



-10/3. This value is changed by in-plane fluid and elastic effects 1 10,1 
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The renormalization group flow of k^s extracted from the one-loop result (||) has no nontrivial fixed point. If this conclusion 
persists to all orders in perturbation theory, it would imply the absence of a smooth phase with long-ranged correlations. The 
smooth appearance of lipid vesicles in the laboratory can then only be explained by their very large persistence length. An 
alternative explanation has recently been proposed in JT^], where it was argued that the neglected higher order terms in the 
energy (jl|) may give rise to a nonperturbative mechanism, by which the crumpled phase can go over into a smooth phase via a 
sequence of two Kosterlitz-Thouless phase transitions or a single first-order one. 

The renormalization of the surface tension has also been investigated by several authors. The results can be summarized by 
the formula 

/icff = /io + i^-ln(AL), (4) 



with the value a' = 1 found in [|^,[P7|] and a' = 3 in [0^]. In Ref. [[L8[], an attempt was made to reconcile the differences. 
An almost planar surface without overhangs was considered in the Monge parametrization. The points on the surface are then 
specified by a vector field r(a;) = {xi,X2, <p{x)), where (f){x) denotes the vertical displacement of the surface with respect to 
a base plane with Cartesian coordinates x = (xi,X2)- For a surface with fixed topology, the Gaussian curvature energy is 
a constant and can be ignored. The remainder of Eq. (|l]) was expanded to fourth order in the displacement field. The relative 
weights of the resulting terms are fixed by their covariant origin. The authors encountered considerable problems in showing that 
this remains true after including the thermal fluctuations. They studied the renormalization of the surface tension by determining 
the coefficient r of the first (constant) term in the expansion of the surface energy. 



1 f (Pk 
^° ^ 2/3 J (2^ ln(/^ofc^ + Kok'^) 



= A'o + T^— In(AL) + ciA2 + 02^^ hi A + cg, (5) 

47rp Kq 

with a constants, and comparing it with the coefficient ^cff of the second term [proportional to 



1 r d^fc 



^j,ok^ + Kok'^ 



= ^^0 + 7^— ln(AL) + C4A2 + c5. (6) 

The above expressions differ from each other by positive powers of the cutoff A. Since the covariance of the theory implies 
T = /icff, two factors were added to the energy, in order to correct Eq. (^. The first one corresponds to the Faddeev-Popov 
determinant associated with fixing the gauge. In the Monge gauge, its contribution is proportional to the cutoff A. The second 
(more ad hoc) factor introduces a nonlinear correction to the integration measure in the partition function. It accounts for the 
difference between an infinitesimal surface element on the membrane and its projection on the reference plane. This second 
correction factor, too, contributes only with positive powers of the cutoff. Added to the first one, it leads to the equality = t, 
which is the main result of Ref. [ jl^ ] . 

Motivated by these problems, and by the renewed interest in the subject, we study the role of thermal fluctuations in a more 
general approach. Employing a derivative expansion [[l9|], we calculate the full effective energy functional produced by Gaussian 
fluctuations for an arbitrary background configuration, maintaining the full nonlinear structure of the energy at all intermediate 
steps. 

The mean curvature in the Monge parametrization reads 

H = d-N^d^,N^, (7) 
where the summation is over the first two components only (/i = 1, 2); N is the unit normal to the surface 

\ (-910,-^20,1), (8) 

y/1 + (90)2 

and the surface elements are 



dS = d^x^l + (90)2, (9) 
so that the first two terms of the energy (|l]) read explicitly 
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1 + (90)2 
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,di,(f)d'^(j) {df,<j>d^(l)df,d^(j>f 
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Physically, /io corresponds to the chemical potential specifying the exchange of molecules between the (incompressible) mem- 
brane and its aqueous environment. 

The main purpose of this note is to show that the ultraviolet divergent parts of the one-loop corrections induced by thermal 
fluctuations are of precisely the same form as in (p^, and in particular, that the three terms in the curvature energy renormalize 
in the same way, resulting in an overall renormalization of kq alone. 

To apply the derivative expansion we write the partition function as a functional integral over the displacement field 



Z 



(11) 



with each field configuration weighted with a Boltzmann factor. Fixing a gauge is generally accompanied by a Faddeev-Popov 
determinant appearing in the measure of the functional integral. Following Ref. we adopt dimensional regularization to 
handle momentum integrals which diverge in the ultraviolet. This is common practice in the technically closely related nonlin- 
ear sigma model. The great advantage of dimensional regularization over regularization with a momentum cutoff A ^ 1/a, is 
that terms diverging with a strictly positive power of A are suppressed. As a result, both the Faddeev-Popov determinant corre- 
sponding to the Monge gauge and the second correction factor introduced in Ref. [jl^, which contain only positive powers of 
the cutoff, are unity in dimensional regularization, and the difficulties addressed in that reference are avoided. Only logarithmic 
divergences show up as poles in e, where e — 2 — D, D being the dimension of the membrane. The connection between the two 
types of regularization is 
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In(AL); 



(12) 



with the linear size L of the membrane representing the relevant long-distance scale. The rationale for using dimensional 
regularization is that contributions to the effective energy with strictly positive powers of the ultraviolet cutoff are connected to 
(5^2) (a; — 0). These highly local terms are uninteresting at large lenght scales [ pTp2| ]. 

In the one-loop approximation, the exponent in ( pi] ) may be expanded up to second order around a background configuration 
$(a;) extremizing Eq. A nontrivial background requires the presence of an extra source term. For brevity, this term will not be 
written down explicitly when setting SEq = 0. The resulting integral is Gaussian and yields an effective energy 



E.si'P] ^ Eom + E^m ^ Eo[^] 
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where the expression in square brackets corresponds to the matrix of second functional derivatives of Eq and the trace Tr stands 
for the trace of this matrix, i.e., the integral J d^x over space, as well as the integral / d2fc/(27r)2 over momentum [p^. 

The one-loop correction Ei [$] to the energy will now be calculated in a derivative expansion for a nearly flat, but arbitrarily 
tilted background configuration. The expansion has the general form 

E,[^] = J d^x [V{Vx) + Z\Vx){d,,V^)^ + Zl^{Vx)d^V,d„V„ + Zl^^p{Vx)d^,V,d„Vp + •••], (14) 

where we introduced the abbreviation V^i, — 9^$, while V, Z^, Z^^, and Z^^^^.^ are functions of to be determined. FoHowing 
Ref. [p^, we set V^i(a;) = + Vf^{x), where V"^ denotes the constant part of V^(x), and expand Eq. ( [l4| ) in powers of Vfj_{x) 
and its derivatives, to obtain 



Ei[Vx+vx] 



(15) 



with space-independent V{V\) and Z(Vx)'s. These functions will now be extracted from the expansion of the Tr In in ( 13 ) up to 
quadratic terms in and d^j^Vu. 

The functional derivatives in (|l3|) are calculated using the Euler-Lagrange formula 
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To keep track of the many terms appearing in the resulting expression we have used 



with F[(f>] = J d^xfi^, d^^, d^d^cj), . . 

the algebraic computer program FORM 

We consider first the renormalization of the surface tension. Since the energy density vl + V^ does not contain derivatives 
of V^, we may set v^ {x) to zero and consider Ei \V\] only, 



l3Ei[Vx] = -iTrln(l + V^) + \lx\n[G-\p)]. 
Here, G~^{p) denotes the inverse propagator: 

G-'ip) = (mo + Atop')p' - (/io + 2KoP^){U-pf + KoitJ-pY 



where is the constant vector 



(17) 



(18) 



(19) 



In dimensional regularization, the first term at the right-hand side of ( |17| ) is zero. To evaluate the remaining Tr In, we apply a 
standard trick and first differentiate (|l7|) with respect to no to obtain 



P 



[U-pf 



— Tr 

2/3 ' L(/^o + A«op2)p2 - (mo + 2kop^){U ■ pf + Ko{U- pY 



where p^ = —id^. Because the integrand contains no space-dependence, the spatial part of the trace in (^ 
A ^ J d^x, and we are left with the momentum integral 



dEijVx 
dfjQ 



A 



d^k 



fc2 - {u-ky 



(no + Kofc2)fc2 _ + 2Kok^){U ■ kf + Ko(f7 • ky 



(20) 

I yields an area factor 
(21) 



Being interested only in the ultraviolet divergent terms, we obtain, in dimensional regularization 



dE^[Vx 



1 1 



d^x^/l + V^ 



After integrating again with respect to fiQ and comparing the result with ( p^ we find (up to an irrelevant additive constant) 



V(14) 



47r/3Ko e 



Vi + (W, 



(22) 



(23) 



where we replaced Vx with the full background field Vx{x), to obtain the first term in (|l^. Note that this one-loop correction is 
precisely of the same form as the surface term contained in the original energy expression (|]). This term can consequently be 
combined with the original one by introducing the renormahzed tension 



fJ-cS = Mo + 



1 Mo 1 

An (3 kq e 



(24) 



This result, corresponding to a' = 1 in (^, is in agreement with Refs. [plJlTl], but disagrees with Refs. [|7|,^ where the value 
a' — S was obtained. To understand the differences, we note that in these last two references, the energy (|l]) with mq = 
was used instead. That is, the renormalization of the surface tension calculated by these authors was generated solely by the 
curvature terms. However, the surface term also contributes. In fact, it generates a contribution with a' = —2, which, together 
with the contribution obtained from the curvature terms, results in the value a' = 1. This is also the value obtained in Ref. JT^], 
as can be seen by disregarding the terms proportional to positive powers of A in Eq. (^ and using the relation (|l2|). There, the 
covariance of the first two terms in the expansion of the surface energy is assured by introducing correction factors proportional 
to positive powers of the cutoff. Our result, based on dimensional regularization where terms with positive powers of the cutoff 
are suppressed, proves the covariance of all terms in the expansion of the surface energy since the full expression has been 
maintained. 

We continue to investigate the renormalization of the bending rigidity. Since the three terms involved contain derivatives of 
the background field V^, we now have to employ the derivative expansion. As a first step, we Taylor expand the logarithm in 
(111) as: 

l3Ei[Vx + vx{x)] - pEi[Vx] = iTrln[l + G{p)K{x,p)] = iTr[G(p)A(p, x)] - \i:r[G{p)K{x,p)G{p)K{x,p)] + • • • , (25) 
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where G{p) is the propagator defined in and A(a;, p) contains the x-dependent terms obtained from functionally differenti- 
ating Eq twice, setting ^^^(x) = V^{x) = + v^{x) and expanding up to second order in and df^v„. 

The first term in (|2q) can be calculated in a similar fashion as V{V\). In the second term, all momentum operators have to be 
moved to the left [|19|], by repeatedly applying the identity 



fix)Pt,g{x) = {pf, + idf,)f{x)g{x), 



(26) 



where f{x) and g{x) are arbitrary functions and the derivative 9^ acts only on the next object to the right, while the derivative 
Pf^ acts on everything to the right. 

The typical momentum integrals showing up at the one-loop order are of the form 



{2n) 



D 



dfc 



infrared 
ultraviolet ' 



(27) 



with m,n > 0. They diverge in the infrared when m + D— l~2n< —1, and in the ultraviolet when m + D — 1 — An> — 1. 
For D = 2 these conditions become to — 2n < — 2, to — 4n > —2, respectively, and the two types of divergences are seen to be 
separated by a wedge of finite integrals in the (m, n)-plane starting at (—2, 0). 

After a tedious and lengthy calculation, involving of the order of 10* terms, done with help of a program written in FORM 
p^], we obtained the divergent terms to second order in derivatives of the field Vfj,: 



(3E,[Vx + vx] - l3E,[Vx] ^ / d'x 



Attkq e 
3 1 
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An ew 
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(1-1-172)1/2 

(i + y2)3/ 



(l-f 1/2)1/2 (1 + ^2)3/2 



'(1 + ^2)3/2 ^ 



y2)5/2 



(l + y2)5/ 




dpVudcrVp 



(28) 



In deriving this expression we also encountered infrared divergences. These are regularized in the same scheme as used to 
regularize the ultraviolet divergences. To distinguish the two we gave epsilon an index ir in case of an infrared divergence. We 
leave the discussion of the infrared divergences to the next paragraph, and first analyze the ultraviolet ones. Comparing ( ^8| ) 
to ( p^ ) with V{V\) given by (p3|), we see that the terms proportional to /io precisely correspond to the first two terms at the 
right-hand side of (ttSh, as it should be. Moreover, we conclude that the Z-functions in (nsl) are given by 



3 1 



Btt/? e (1 y2)i/; 



3 1 V,V. „3 (y. 
) Zuyap{V\) 



Att/S e (1 + y2)3/2 



3 1 F^KK^p 
'87r/3e (H-y2)5/2- 



(29) 



By replacing the constant V\ with the full background field (x) — dp^{x), we obtain for the divergent parts of the expansion 
( |l4| ) the explicit form 



(C)2$)2 
1 + (9$)' 



- 2 



[1 + (9$)2 



[1 + (a$) 
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(30) 



We see that the thermally generated terms at the one-loop level are precisely of the same form as those present in the original 
energy expression ([lo|). In addition, the relative weights of the curvature terms produced by the fluctuations are the same as 
those found there. They can therefore be combined with the original terms by introducing the renormalized rigidity 



3 1 

We' 



(31) 



whose value is in agreement with [|7[-|11| . 

As seen in (^8|), the one-loop corrections seem to have introduced infrared divergences in the theory. A closer inspection 
reveals that the infrared-divergent contributions all stem from the surface energy term in (p^, so that it suffices to analyze the 
one-loop corrections to the truncated energy 



(32) 
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Infrared divergences in this model have previously been studied in where they were shown to disappear for an infinitely 
small dimension D of the membrane to all orders in D. In our calculation the problem arises for D = 2 — e. When calculating 
the effective action, we expand around the background field $ extremizing E'q, i.e.. 



5E', 



0, (33) 



which reads explicitly 



[i + (a$)2]i/2 [i + (a$)2]3/2 • ^ ' 

The presence of the implicitly assumed sources turns this equation in a nontrivial one. Rewriting d^^{x) = + v^{x), 
expanding to linear order in v^, and substituting the resulting expression in (p^, we see the infrared divergences to vanish for a 
two-dimensional membrane. 

In conclusion, we have demonstrated that all logarithmically divergent one-loop corrections induced by thermal fluctuations 
are precisely of the same form as in the original energy (|lO|), so that they can be removed by a renormalization of the surface 
tension and bending rigidity. 
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